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Abstract 


Semantic word embeddings represent the meaning of a word via a vector, and are created by diverse 
methods. Many use nonlinear operations on co-occurrence statistics, and have hand-tuned hyperparam¬ 
eters and reweighting methods. 

This paper proposes a new generative model, a dynamic version of the log-linear topic model of|Mnih| 


and Hinton (20071. The methodological novelty is to use the prior to compute closed form expressions 
for word statistics. This provides a theoretical justification for nonlinear models like PMI, word2vec, 
and GloVe, as well as some hyperparameter choices. It also helps explain why low-dimensional semantic 
embeddings contain linear algebraic structure that allows solution of word analogies, as shown by |Mikolov| 
et al.| ( |2013a[ ) and many subsequent papers. 

Experimental support is provided for the generative model assumptions, the most important of which 
is that latent word vectors are fairly uniformly dispersed in space. 


1 Introduction 


Vector representations of words (word embeddings) try to capture relationships between words as distance or 
angle, and have many applications in computational linguistics and machine learning. They are constructed 
by various models whose unifying philosophy is that the meaning of a word is defined by “the company it 
keeps” (Firth 1957| , namely, co-occurrence statistics. The simplest methods use word vectors that explicitly 
represent co-occurrence statistics. Reweighting heuristics are known to improve these methods, as is dimen¬ 


sion reduction (Deerwester et al. 1990). Some reweighting methods are nonlinear, which include taking the 


square root of co-occurrence counts (Rohde et al. 2006), or the logarithm, or the related Pointwise Mutual 


Information (PMI) (Church and Hanks 1990). These are collectively referred to as Vector Space Models, 
surveyed in ([Turney and Pantelf 2010). 


Neural network language models (Rumelhart et al. 1986 1988 Bengio et al. 2006 Collobert and Weston 


2008a) propose another way to construct embeddings: the word vector is simply the neural network’s internal 


representation for the word. This method is nonlinear and nonconvex. It was popularized via word2vec, a 
family of energy-based models in (Mikolov et al. 2013b|c ), followed by a matrix factorization approach called 
GloVe (Pennington et al. 2014). The first paper also showed how to solve analogies using linear algebra on 


word embeddings. Experiments and theory were used to suggest that these newer methods are related to 


the older PMI-based models, but with new hyperparameters and/or term reweighting methods (Levy and 
Goldberg 2014b). 


But note that even the old PMI method is a bit mysterious. The simplest version considers a symmetric 
matrix with each row/column indexed by a word. The entry for {w, w') is PMI(t(;, w') — log “tAwA, where 


p{w^ w') is the empirical probability of words w, w' appearing within a window of certain size in the corpus, 
and p{w) is the marginal probability of w. (More complicated models could use asymmetric matrices with 
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columns corresponding to context words or phrases, and also involve tensorization.) Then word vectors are 
obtained by low-rank SVD on this matrix, or a related matrix with term reweightings. In particular, the 
PMI matrix is found to be closely approximated by a low rank matrix: there exist word vectors in say 300 
dimensions, which is much smaller than the number of words in the dictionary, such that 

~ ( 1 . 1 ) 


where « should be interpreted loosely. 

There appears to be no theoretical explanation for this empirical finding about the approximate low 
rank of the PMI matrix. The current paper addresses this. Specifically, we propose a probabilistic model 


of text generation that augments the log-linear topic model of Mnih and Hinton (20071 with dynamics, in 


the form of a random walk over a latent discourse space. The chief methodological contribution is using the 


model priors to analytically derive a closed-form expression that directly explains (1.1); see Theorem 2.2 


in Section Section builds on this insight to give a rigorous justification for models such as word2vec 
and GloVe, including the hyperparameter choices for the latter. The insight also leads to a mathematical 
explanation for why these word embeddings allow analogies to be solved using linear algebra; see Section]^ 
Section shows good empirical fit to this model’s assumtions and predictions, including the surprising one 
that word vectors are pretty uniformly distributed (isotropic) in space. 


1.1 Related work 

Latent variable probabilistic models of language have been used for word embeddings before, including 


Latent Dirichlet Allocation (LDA) and its more complicated variants (see the survey (Blei 2012)), and some 
neurally inspired nonlinear models ( |Mnih and Hinton 2007 Maas et al. 2011). In fact, LDA evolved out of 
efforts in the 1990s to provide a generative model that “explains” the success of older vector space methods 


like Latent Semantic Indexing (Papadimitriou et al. 1998 Hofmann 1999). However, none of these earlier 


generative models has been linked to PMI models. 


Levy and Goldberg (2014b) tried to relate word2vec to PMI models. They showed that if there were no 


dimension constraint in word2vec, specifically, the “skip-gram with negative sampling (SGNS)” version of the 
model, then its solutions would satisfy (1.11, provided the right hand side were replaced by PMI(w, w') — /3 
for some scalar (3. However, skip-gram is a discriminative model (due to the use of negative sampling), not 
generative. Furthermore, their argument only applies to very high-dimensional word embeddings, and thus 
does not address low-dimensional embeddings, which have superior quality in applications. 


Hashimoto et al. (2016) focuses on issues similar to our paper. They model text generation as a random 


walk on words, which are assumed to be embedded as vectors in a geometric space. Given that the last word 
produced was w, the probability that the next word is w' is assumed to be given by h{\vuj — Uu,'P) for a 


suitable function h, and this model leads to an explanation of (1.1). By contrast our random walk involves 


a latent discourse vector, which has a clearer semantic interpretation and has proven useful in subsequent 
work, e.g. understanding structure of word embeddings for polysemous words Arora et al.| ( |2016[ ). Also our 
work clarifies some weighting and bias terms in the training objectives of previous methods (Section]^ and 
also the phenomenon discussed in the next paragraph. 

Researchers have tried to understand why vectors obtained from the highly nonlinear word2vec models 


exhibit linear structures (Levy and Goldberg 2014a Pennington et al., 2014). Specifically, for analogies like 


^‘man:woman::king:13 ” queen happens to be the word whose vector Vqueen is the most similar to the vector 
Vking — Vman + Vwoman- This suggests that simple Semantic relationships, such as masculine vs feminine 
tested in the above example, correspond approximately to a single direction in space, a phenomenon we will 
henceforth refer to as RELATI0NS=LINES. 

Section surveys earlier attempts to explain this phenomenon and their shortcoming, namely, that they 
ignore the large approximation error in relationships like (1.1). This error appears larger than the difference 


between the best solution and the second best (incorrect) solution in analogy solving, so that this error could 
in principle lead to a complete failure in analogy solving. In our explanation, the low dimensionality of the 
word vectors plays a key role. This can also be seen as a theoretical explanation of the old observation that 
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dimension reduction improves the quality of word embeddings for various tasks. The intuitive explanation 
often given —that smaller models generalize better—turns out to be fallacious, since the training method 
for creating embeddings makes no reference to analogy solving. Thus there is no a priori reason why low¬ 
dimensional model parameters (i.e., lower model capacity) should lead to better performance in analogy 
solving, just as there is no reason they are better at some other unrelated task like predicting the weather. 


1.2 Benefits of generative approaches 

In addition to giving some form of “unification” of existing methods, our generative model also brings 
more intepretability to word embeddings beyond traditional cosine similarity and even analogy solving. For 
example, it led to an understanding of how the different senses of a polysemous word (e.g., bank) reside in 
linear superposition within the word embedding (Arora et al., 2016). Such insight into embeddings may 


prove useful in the numerous settings in NLP and neuroscience where they are used. 

Another new explanatory feature of our model is that low dimensionality of word embeddings plays a key 
theoretical role —unlike in previous papers where the model is agnostic about the dimension of the embed¬ 


dings, and the superiority of low-dimensional embeddings is an empirical hnding (starting with Deerwester 
et al.| (1990)). Specifically, our theoretical analysis makes the key assumption that the set of all word vectors 
(which are latent variables of the generative model) are spatially isotropic, which means that they have no 
preferred direction in space. Having n vectors be isotropic in d dimensions requires d n. This isotropy is 


needed in the calculations (i.e., multidimensional integral) that yield (1.1). It also holds empirically for our 
word vectors, as shown in Section 

The isotropy of low-dimensional word vectors also plays a key role in our explanation of the RELA- 
TiONS=LiNES phenomenon (Section]^. The isotropy has a “purihcation” effect that mitigates the effect of 
the (rather large) approximation error in the PMI models. 


2 Generative model and its properties 

The model treats corpus generation as a dynamic process, where the t-th word is produced at step t. The 
process is driven by the random walk of a discourse vector Ct € 3?'^. Its coordinates represent what is being 
talked aboutj^ Each word has a (time-invariant) latent vector Vw € 3?^ that captures its correlations with 
the discourse vector. We model this bias with a log-linear word production model: 


Pr[ri; emitted at time t | c*] oc exp((ct,Vw})- 


( 2 . 1 ) 


The discourse vector ct does a slow random walk (meaning that ct+i is obtained from ct by adding a small 
random displacement vector), so that nearby words are generated under similar discourses. We are interested 
in the probabilities that word pairs co-occur near each other, so occasional big jumps in the random walk 
are allowed because they have negligible effect on these probabilities. 

A similar log-linear model appears in Mnih and Hinton ( 2007| but without the random walk. The linear 
chain CRF of Collobert and Weston (2008b) is more general. The dynamic topic model of Blei and Lafferty 


(2006) utilizes topic dynamics, but with a linear word production model. Belanger and Kakade (2015) have 


proposed a dynamic model for text using Kalman Filters, where the sequence of words is generated from 
Gaussian linear dynamical systems, rather than the log-linear model in our case. 


The novelty here over such past works is a theoretical analysis in the method-of-moments tradition (Hsu 


et al. 

2012 

Cohen et al. 

2012 


Assuming a prior on the random walk we analytically integrate out the 
hidden random variables and compute a simple closed form expression that approximately connects the 


model parameters to the observable joint probabilities (see Theorem 2.2). This is reminiscent of analysis of 


similar random walk models in finance (Black and Scholes 1973). 


^This is a different interpretation of the term “discourse” compared to some other settings in computational linguistics. 
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Model details. Let n denote the number of words and d denote the dimension of the discourse space, 
where 1 < d < n. Inspecting (2.1) suggests word vectors need to have varying lengths, to fit the empirical 


finding that word probabilities satisfy a power law. Furthermore, we will assume that in the bulk, the word 
vectors are distributed uniformly in space, earlier referred to as isotropy. This can be quantified as a prior 
in the Bayesian tradition. More precisely, the ensemble of word vectors consists of i.i.d draws generated by 
V = s ■ V, where v is from the spherical Gaussian distribution, and s is a scalar random variable. We assume 
s is a random scalar with expectation r = 0(1) and s is always upper bounded by k, which is another 
constant. Here r governs the expected magnitude of {v,Ct), and it is particularly important to choose it to 
be 0(1) so that the distribution Pr[i(;|ct] cx exp((?;„,, ct)) is interesting!^ Moreover, the dynamic range of 
word probabilities will roughly equal exp(K^), so one should think of k as an absolute constant like 5. These 
details about s are important for realistic modeling but not too important in our analysis. (Furthermore, 
readers uncomfortable with this simplistic Bayesian prior should look at Section 2.2 below.) 

Finally, we clarify the nature of the random walk. We assume that the stationary distribution of the 
random walk is uniform over the unit sphere, denoted by C. The transition kernel of the random walk can 
be in any form so long as at each step the movement of the discourse vector is at most in li normj^ 

This is still fast enough to let the walk mix quickly in the space. 

The following lemma (whose proof appears in the appendix) is central to the analysis. It says that 
under the Bayesian prior, the partition function c)), which is the implied normalization 

in equation (2.1), is close to some constant Z for most of the discourses c. This can be seen as a plausible 


theoretical explanation of a phenomenon called self-normalization in log-linear models: ignoring the partition 
function or treating it as a constant (which greatly simplifies training) is known to often give good results. 
This has also been studied in (Andreas and Klein[ 2014). 


Lemma 2.1 (Concentration of partition functions). If the word vectors satisfy the Bayesian prior described 
in the model details, then 

Pr [(1 - e,)Z <Z,<{1 + e,)Z] >1-6, (2.2) 


c~C 


for €z = 0{l/ydn), and 6 = exp(—H(log^ n)). 

The concentration of the partition functions then leads to our main theorem (the proof is in the appendix). 
The theorem gives simple closed form approximations for p{w), the probability of word w in the corpus, and 
p(w,w'), the probability that two words w,w' occur next to each other. The theorem states the result for 
the window size q = 2, but the same analysis works for pairs that appear in a small window, say of size 10, 


as stated in Corollary 2.3 Recall that PMI(i(;,r(;') = \og[p{w, w')/{p{w)p{w'))]. 

Theorem 2.2. Suppose the word vectors satisfy the inequality (2.2), and window size q = 2. Then, 


\ogp{w, w') = - 2 log Z ± e. 


logp(r/;) = 


2d 


— log Z ± e. 


(2.3) 

(2.4) 


for e = 0{ez) -b 0{l/d) -I- 0(e2)- Jointly these imply: 


PMI{w,w') = ±Q(e). (2.5) 

Remarks 1. Since the word vectors have (2 norm of the order of Vd, for two typical word vectors Vw,Vw', 
ll^u, + is of the order of Q[d). Therefore the noise level e is very small compared to the leading 

larger r will make Pr[R;|ct] too peaked and a smaller one will make it too uniform. 

^ More precisely, the proof extends to any symmetric product stationary distribution C with sub-Gaussian coordinate 
satisfying Ec [l|c||^] = 1 , and the steps are such that for all Cf, Ep(ci_^i|c£)— ct||)] <1 + 62 for some small £ 2 - 
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term i^\\vw + Wiu'lli- however, the noise level 0{e) could be comparable to the leading term, and 

empirically we also find higher error here. 


Remarks 2. Variants of the expression for joint probability in (2.3) had been hypothesized based upon 


empirical evidence in Mikolov et al. (2013b) and also Globerson et al. (2007), and Maron et al. (2010) . 


Remarks 3. Theorem |2.2| directly leads to the extension to a general window size q as follows: 

Corollary 2.3. Let pq{w,w') be the co-occurrence probability in windows of size q, and PMIq(w,w') be the 
corresponding PMI value. Then 

logpq{w,w') = - 21 ogZ + 7 ±e, 

PMIq {w,w') = + 7 ±C>(e). 


where 7 = log ^ 


9(g-l) 

2 


It is quite easy to see that Theorem |2.2| implies the Corollary |2.3[ as when the window size is q the 


pair w^w' could appear in any of (|) positions within the window, and the joint probability of w,w' is 
roughly the same for any positions because the discourse vector changes slowly. (Of course, the error 
term gets worse as we consider larger window sizes, although for any constant size, the statement of the 


theorem is correct.) This is also consistent with the shift /? for fitting PMI in (Levy and Goldberg 2014b), 


which showed that without dimension constraints, the solution to skip-gram with negative sampling satisfies 
PMI {w,w') — /? = {vw^Vw') for a constant P that is related to the negative sampling in the optimization. 
Our result justifies via a generative model why this should be satisfied even for low dimensional word vectors. 


2.1 Proof sketches 

Here we provide the proof sketches, while the complete proof can be found in the appendix. 


Proof sketch of Theorem 2.2 Let w and w' be two arbitrary words. Let c and c' denote two consecutive 
context vectors, where c ^ C and c'|c is dehned by the Markov kernel p{c' \ c). 

We start by using the law of total expectation, integrating out the hidden variables c and cb 

p{w,w') = E [Pr[r(;,w'lc,c']] 

c,c' 

= E [p{w\c)p{w'\c')] 

c,c' 

'exp((uu,,c)) exp((uu,/,c')) 


= E 


( 2 . 6 ) 


An expectation like (2.6) would normally be difficult to analyze because of the partition functions. 


However, we can assume the inequality (2.2), that is, the partition function typically does not vary much for 


most of context vectors c. Let F be the event that both c and c' are within (l±e^)Z. Then by (2.2) and the 
union bound, event F happens with probability at least 1 — 2exp(—H(log^ n)). We will split the right-hand 
side (RHS) of (2.6) into the parts according to whether F happens or not. 


RHS of (2.6) = E 


exp((ti„,,c)) exp((u^,/,c')) ^ 

Zr. Zr' 


E 


exp((i 


Ti 

, c)) exp{{v^ 


c')). 


(2.7) 


T2 
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where F denotes the complement of event F and Ip and 1^ denote indicator functions for F and F, 
respectively. When F happens, we can replace Zc by Z with a 1 ± Cz factor loss: The first term of the RHS 
of (2.7) equals to 


Ti = ^ E [exp((i;u,, c)) exp(c'))I f] 

c.c' 


( 2 . 8 ) 


On the other hand, we can use E[1 f] = Pr[.F] < exp(—r2(log^ n)) to show that the second term of RHS 
of (2.7) is negligible, 

|T 2 I = exp(-H(log^-® n)). (2.9) 

This claim must be handled somewhat carefully since the RHS does not depend on d at all. Briefly, the 
reason this holds is as follows: in the regime when d is small {Vd = o(log^n)), any word vector and 
discourse c satisfies that exp((r;u,,c)) < exp(||z)u,||) = exp(0(-\/d)), and since E[1 f] = exp(—H(log^ n)), the 
claim follows directly; In the regime when d is large (\/d = r2(log^ n)), we can use concentration inequalities 
to show that except with a small probability exp(—H(d)) = exp(—H(log^n)), a uniform sample from the 
sphere behaves equivalently to sampling all of the coordinates from a standard Gaussian distribution with 
mean 0 and variance in which case the claim is not too difficult to show using Gaussian tail bounds. 

Therefore it suffices to only consider (2.8). Our model assumptions state that c and c' cannot be too 
different. We leverage that by rewriting (2.8) a little, and get that it equals 




Z^ 

1 ± 0(ez) 


exp((nu,,c)) E [exp((•(;„/, c'))] 

c'\c 


E[exp(('Uu,,c))H(c)] 


( 2 . 10 ) 


where A{c) := E [exp((z;u,/, c'))]. We claim that A{c) = (1 ± 0(e2)) exp((ti„/, c)). Doing some algebraic 

c'|c 

manipulations, 

A{c) = exp((vu,',c)) E [exp((z;u,/,c' - c))]. 

c'\c 

Furthermore, by our model assumptions, jjc— c'jj < til'/d- So 

{Vn,,c- c') < ||riu,||||c- c'll = 0 ( 62 ) 


and thus A{c) = (1 ± 0(e2)) exp(c)). Plugging the simplification of A{c) to (2.10), 

l±0(e.) 


Ti = 


Z2 


E[exp(('Uu, + •(;„/, c))]. 


( 2 . 11 ) 


Since c has uniform distribution over the sphere, the random variable {vw +Vu!',c) has distribution pretty 
similar to Gaussian distribution Af{0, ||wu, + Vuj'W^/d), especially when d is relatively large. Observe that 
E[exp(X)] has a closed form for Gaussian random variable X ^ A/'(0, cr^), 

f I 

E[exp(X)] = J —^ exp(a;)dx 

= exp(cr^/2). (2.12) 

Bounding the difference between {v^ + v^', c) from Gaussian random variable, we can show that for 

e = 0(l/d). 


E[exp((nu, + Vyj',c))] = (1 ± e) exp 


\Vw + V-u 

2d 


(2.13) 


Therefore, the series of simplification/approximation above (concretely, combining equations (2.6), (2.7), (2.9), (2.11), 
and ( |2.13 )) lead to the desired bound on \ogp{w,w') for the case when the window size q = 2. The bound 
on logp(?c) can be shown similarly. 
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Proof sketch of Lemma |2.1| Note that for fixed c, when word vectors have Gaussian priors assumed as 
in our model, exp((uu,, c)) is a sum of independent random variables. 

We hrst claim that using proper concentration of measure tools, it can be shown that the variance of 
are relatively small compared to its mean thus Z^ concentrates around its mean. Note this 

is quite non-trivial: the random variable exp((r!u,, c)) is neither bounded nor subgaussian/sub-exponential, 
since the tail is approximately inverse poly-logarithmic instead of inverse exponential. In fact, the same 
concentration phenomenon does not happen for w. The occurrence probability of word w is not necessarily 
concentrated because the £2 norm of can vary a lot in our model, which allows the frequency of the words 
to have a large dynamic range. 

So now it suffices to show that \Z^ for different c are close to each other. Using the fact that the 
word vector directions have a Gaussian distribution, E,;„ [Zc] turns out to only depend on the norm of c 
(which is equal to 1). More precisely, 

E[^c] = /(||c||^) = /(l) (2.14) 

where / is defined as f(a) = nEs[exp(s^a/2)] and s has the same distribution as the norms of the word 
vectors. We sketch the proof of this. In our model, = Sm-Vm, where is a Gaussian vector with identity 
covariance I. Then 


E [Z^j = n E [exp(('Uu,, c))] 


= n E 


E 

Vw I Si 


[exp((u^,,c)) I 


where the second line is just an application of the law of total expectation, if we pick the norm of the 
(random) vector first, followed by its direction. Gonditioned on s^,, (uu,,c) is a Gaussian random variable 
| 2 S^, and therefore using similar calculation as in (2.12|), we have 


E 

™|s. 


[exp((uu,,c)) I s„] = exp(s^||c|| 2 / 2 ). 


Hence, Ev„lZc] = nEs[exp(s^||c||2/2)] as needed. 


Proof of Theorem 4.1 The proof uses the standard analysis of linear regression. Let V = be the 


SVD of V and let cti, ..., cr^ be the left singular values of V (the diagonal entries of E). For notational ease 
we omit the subscripts in ^ and since they are not relevant for this proof. Since and thus 

( = V^C = QE-ip^C', we have 


We claim 


||Cl2<a,-i||P^C'l|2. 



(2.15) 


(2.16) 


Indeed, = 0{nd), since the average squared norm of a word vector is d. The claim then follows 

from the first assumption. Furthermore, by the second assumption, ||P^C1loo < ;^IIC^Il2i so 


p^cr2<-^\\c 


'l|2 

2- 


Plugging (2.161 and (2.17) into (2.15), we get 


IICI 2 < 



C2\/d 

y/c{n 


IIC'I|2 


(2.17) 


as desired. The last statement follows because the norm of the signal, which is dlog(^fl) originally and is 
d\og{vji) = Va — Vb after dimension reduction, also gets reduced by a factor of -^/n. 
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2.2 Weakening the model assumptions 

For readers uncomfortable with Bayesian priors, we can replace our assumptions with concrete properties of 
word vectors that are empirically verifiable (Section |5.1[ ) for our final word vectors, and in fact also for word 
vectors computed using other recent methods. 

The word meanings are assumed to be represented by some “ground truth” vectors, which the experi¬ 
menter is trying to recover. These ground truth vectors are assumed to be spatially isotropic in the bulk, 
in the following two specific ways: (i) For almost all unit vectors c the sum c)) is close to a 

constant Z; (ii) Singular values of the matrix of word vectors satisfy properties similar to those of random 
matrices, as formalized in the paragraph before Theorem |4.1| Our Bayesian prior on the word vectors hap¬ 
pens to imply that these two conditions hold with high probability. But the conditions may hold even if the 
prior doesn’t hold. Furthermore, they are compatible with all sorts of local structure among word vectors 
such as existence of clusterings, which would be absent in truly random vectors drawn from our prior. 


3 Training objective and relationship to other models 


To get a training objective out of Theorem |2.2[ we reason as follows. Let X^i^w' be the number of times words 
w and w' co-occur within the same window in the corpus. The probability p{w,w') of such a co-occurrence 
at any particular time is given by (2.3). Successive samples from a random walk are not independent. 


But if the random walk mixes fairly quickly (the mixing time is related to the logarithm of the vocabulary 
size), then the distribution of is very close to a multinomial distribution Mul(L, w')}), where 

^ ~ w' ^w,w' is the total number of word pairs. 

Assuming this approximation, we show below that the maximum likelihood values for the word vectors 
correspond to the following optimization. 


min y X. 

W^W 


w,w' 


^log( ) 


-C] 


As is usual, empirical performance is improved by weighting down very frequent word pairs, possibly 
because very frequent words such as “the” do not fit our model. This is done by replacing the weighting 
Xw,w' by its truncation min{V.u,^u,',Xi„ax} where X^ax is a constant such as 100. We call this objective with 
the truncated weights SN (Squared Norm). 

We now give its derivation. Maximizing the likelihood of is equivalent to maximizing 


^ = log ( p{w,w'y 


Denote the logarithm of the ratio between the expected count and the empirical count as 


= log 


Lp{w, w') 


(3.1) 


Then with some calculation, we obtain the following where c is independent of the empirical observations 


X„ 


o' s. 


E 


Y A 

■^w.w' ^w.w' 


{w^w') 
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(3.2) 




On the other hand, using ~ 1 + x + a:^/2 when x is smallj^ we have 


L - ^ ^ Lpyjyj! - ^ ^ 

{w,w') {w,w') 


A„ 


{w,w') \ / 


Note that L = Y.{w,w') so 


A , ~ W Y , a2 

■^w^w' ^w,w' ~ f-j / ^ -^w^w' ‘-^w^w' ’ 


{w,w') {w^w') 


Plugging this into (3.2) leads to 


{w,w') 


(3.3) 


So maximizing the likelihood is approximately equivalent to minimizing the right hand side, which (by 
examining (3.1)) leads to our objective. 


Objective for training with PMI. A similar objective PMI can be obtained from (2.5), by computing 
an approximate MLE, using the fact that the error between the empirical and true value of PMI(t(;,w') is 
driven by the smaller term p{w,w'), and not the larger terms p(w),p(w'). 


min 

M,c 


Xu,,w' (PMl{w,w') - {v^, 

w,w' 


Vw')) 


2 


This is of course very analogous to classical VSM methods, with a novel reweighting method. 

Fitting to either of the objectives involves solving a version of Weighted SVD which is NP-hard, 
empirically seems solvable in our setting via AdaGrad (]Duchi et al. 2011). 


but 


Connection to GloVe. 


Compare SN with the objective used by GloVe ({Pennington et ah 2014): 


^ ^ ./( ^w' G) 

w^w' 


with f{Xu]^w') = niinjAiyy, 100}. Their weighting methods and the need for bias terms Siu,Su)',G were 
derived by trial and error; here they are all predicted and given meanings due to Theorem |2.2[ specifically 


Connection to word2vec(CBOW). The GBOW model in word2vec posits that the probability of a 
word Wk+i as a function of the previous k words wi,W 2 , ■ ■ ■ ^Wk- 


1 ^ 

p (wk+i\ cx exp((z;„^_^i, - 


^This Taylor series approximation has an error of the order of x^, but ignoring it can be theoretically justified as follows. 
For a large , its value approaches its expectation and thus the corresponding yj/ is close to 0 and thus ignoring 


A3 , is well justified. The terms where A^ yj/ is significant correspond to X^j ^/’s that are small. But empirically, X. 


w,w 


obey a power law distribution (see, e.g. [Pennington et al.| < [2014[ |) using which it can be shown that these terms contribute a 
small f raction of the final objective l |3.3| ). So we can safely ignore the errors. Full details appear in the ArXiv version of this 
paper (Arora et al. 2015 b 
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This expression seems mysterious since it depends upon the average word vector for the previous k words. 
We show it can be theoretically justified. Assume a simplified version of our model, where a small window 
of k words is generated as follows: sample c ^ C, where C is a uniformly random unit vector, then sample 
{wi,W 2 ,..., Wk) ^ Vnj., c))IZc. Furthermore, assume Z^ = Z for any c. 


Lemma 3.1. In the simplified version of our model, the Maximum-a-Posteriori (MAP) estimate of c given 
{wi,W 2 ,...,Wk) is II 

Proof. The c maximizing p {c\wi,W 2 , ■ ■ ■ ,Wk) is the maximizer of p(c)p {wi,W 2 , ■ ■ ■ ,Wk\c). Since p{c) = 
p(c') for any c,c', and we have p {wi,W 2 , ■ ■ ■ ,Wk\c) = exp((^j , c))/^, the maximizer is clearly c = 


Thus using the MAP estimate of Ct gives essentially the same expression as CBOW apart from the 
rescaling, which is often omitted due to computational efficiency in empirical works. 


4 Explaining RELATIONS=LINES 


As mentioned, word analogies like “a:6::c:??” can be solved via a linear algebraic expression: 

argmin H'Ca - Vb - Vc Vd\\l, 


(4.1) 


where vectors have been normalized such that ||r ’£;||2 = 1- This suggests that the semantic relationships 
being tested in the analogy are characterized by a straight line [preferred to earlier as relations=lines. 

Using our model we will show the following for low-dimensional embeddings: for each such relation R 
there is a direction pn in space such that for any word pair a, b satisfying the relation, Va — Vb is like p-n 
plus some noise vector. This happens for relations satisfying a certain condition described below. Empirical 
results supporting this theory appear in Section where this linear structure is further leveraged to slightly 
improve analogy solving. 

A side product of our argument will be a mathematical explanation of the empirically well-established su¬ 


periority of low-dimensional word embeddings over high-dimensional ones in this setting (Levy and Goldberg 


2014a). As mentioned earlier, the usual explanation that smaller models generalize better is fallacious. 


We first sketch what was missing in prior attempts to prove versions of relations=lines from first 
principles. The basic issue is approximation error: the difference between the best solution and the 2nd best 


solution to (4.1) is typically small, whereas the approximation error in the objective in the low-dimensional 


solutions is larger. For instance, if one uses our PMI objective, then the weighted average of the termwise 


error in (2.5) is 17%, and the expression in (4.1) above contains six inner products. Thus in principle the 


approximation error could lead to a failure of the method and the emergence of linear relationship, but it 
does not. 


Prior explanations. Pennington et al. (2014) try to propose a model where such linear relationships 


should occur by design. They posit that queen is a solution to the analogy ^^man:woman::king:ll” because 


Pix I king) p{x \ man) 


p{x I queen) p(x \ woman) 


(4.2) 


where p{x \ king) denotes the conditional probability of seeing word x in a small window of text around 
king. Relationship (14.21) is intuitive since both sides will be « 1 for gender-neutral x like “walks’^ or “food !^, 


®Note that this interpretation has been disputed; e.g., it is argued in |Levy and Goldberg] ( |2014a| that | |4.l| can be understood 
using only the classical connection between inner product and word similarity, using which the objective | |4.1| l is slightly 
improved to a different objective called 3COSMUL. However, this “explanation” is still dogged by the issue of large termwise 
error pinpointed here, since inner product is only a rough approximation to word similarity. Furthermore, the experiments in 
Section clearly support the RELATIONS—LINES interpretation. 
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will be > 1 when x is like “he, Henry" and will be < 1 when x is like “dress, she, Elizabeth." This was 


also observed by Levy and Goldberg (2014a|. Given (4.2), they then posit that the correct model describing 
word embeddings in terms of word occurrences must be a homomorphism from (3?^, +) to (3?“'", x), so vector 
differences map to ratios of probabilities. This leads to the expression 

Pw,w' — {’^w: T by,', 

and their method is a (weighted) least squares ht for this expression. One shortcoming of this argument 
is that the homomorphism assumption assumes the linear relationships instead of explaining them from a 
more basic principle. More importantly, the empirical fit to the homomorphism has nontrivial approximation 
error, high enough that it does not imply the desired strong linear relationships. 


Levy and Goldberg (2014b) show that empirically, skip-gram vectors satisfy 

{vw,Vw') « PMI(u;,w') 


(4.3) 


up to some shift. They also give an argument suggesting this relationship must be present if the solution 
is allowed to be very high-dimensional. Unfortunately, that argument does not extend to low-dimensional 
embeddings. Even if it did, the issue of termwise approximation error remains. 


Our explanation. The current paper has introduced a generative model to theoretically explain the 


emergence of relationship (4.3). However, as noted after Theorem 2.2 the issue of high approximation error 


does not go away either in theory or in the empirical fit. We now show that the isotropy of word vectors 


(assumed in the theoretical model and verified empirically) implies that even a weak version of (4.3) is enough 
to imply the emergence of the observed linear relationships in low-dimensional embeddings. 

This argument will assume the analogy in question involves a relation that obeys Pennington et al.’s 


suggestion in (4.2). Namely, for such a relation R there exists function i^_r(-) depending only upon R such 


that for any a, b satisfying R there is a noise function fa.b.Ri') for which: 


p{x I q) 

P{X I b) 


= vrU) ■ L.bMx) 


(4.4) 


For different words x there is huge variation in (4.4), so the multiplicative noise may be large. 


Our goal is to show that the low-dimensional word embeddings have the property that there is a vector 
Pr such that for every pair of words a, b in that relation, Va — Vb = pr -f noise vector, where the noise vector 
is small. 


Taking logarithms of (4.4) results in: 


log 


Pjx I a) 
P{X I b) 


= fog(^'ii(x)) +Ca.h,fl(x) 


Theorem 


2.2 


implies that the left-hand side simplifies to log 




2.2 


e captures the small approximation errors induced by the inexactness of Theorem 
noise! Denoting by V the n x d matrix whose rows are the vectors, we rewrite (4.5) as: 


(4.5) 

- Vb) + ia.bix) where 
This adds yet more 


V{Va - Vb) = dlogivR) + Ca,b,R 


(4.6) 


where log(z//j) in the element-wise log of vector ixr and Cab r~ d{fa,b,R — ^a,b,R) is the noise. 

In essence, (4.6) shows that Va — Ub is a solution to a linear regression in d variables and m constraints. 


with Ca b R bsing ths “noise.” The design matrix in the regression is V, the matrix of all word vectors, which 
in our model (as well as empirically) satisfies an isotropy condition. This makes it random-like, and thus 
solving the regression by left-multiplying by VC the pseudo-inverse of V, ought to “denoise” effectively. We 
now show that it does. 

Our model assumed the set of all word vectors satisfies bulk properties similar to a set of Gaussian vectors. 
The next theorem will only need the following weaker properties. (1) The smallest non-zero singular value 
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(a) SN (b) GloVe (c) CBOW (d) skip-gram 


Figure 1: The partition function The figure shows the histogram of for 1000 random vectors c of 
appropriate norm, as dehned in the text. The cc-axis is normalized by the mean of the values. The values 
for different c concentrate around the mean, mostly in [0.9,1.1]. This concentration phenomenon is predicted 
by our analysis. 

of V is larger than some constant Ci times the quadratic mean of the singular values, namely, ||y||ir/-\/d. 
Empirically we find ci « 1/3 holds; see Section]^ (2) The left singular vectors behave like random vectors 
with respect to Cab w namely, have inner product at most C2\\Cab rW/V^ with Cab some constant C 2 - 

(3) The max norm of a row in V is 0{Cd). The proof is included in the appendix. 

Theorem 4.1 (Noise reduction). Under the conditions of the previous paragraph, the noise in the dimension- 
reduced semantic vector space satisfies 


\\CARh<\\Ca,b.Rh^- 

As a corollary, the relative error in the dimension-reduced space is a factor of \fd[n smaller. 


5 Experimental verification 

In this section, we provide experiments empirically supporting our generative model. 

Corpus. All word embedding vectors are trained on the English Wikipedia (March 2015 dump). It is 
pre-processed by standard approach (removing non-textual elements, sentence splitting, and tokenization), 
leaving about 3 billion tokens. Words that appeared less than 1000 times in the corpus are ignored, resulting 
in a vocabulary of 68,430. The co-occurrence is then computed using windows of 10 tokens to each side of 
the focus word. 


Training method. Our embedding vectors are trained by optimizing the SN objective using AdaGrad (Due hi 


et al. 20111 with initial learning rate of 0.05 and 100 iterations. The PMI objective derived from (2.5) was 


also used. SN has average (weighted) term-wise error of 5%, and PMI has 17%. We observed that SN 
vectors typically fit the model better and have better performance, which can be explained by larger errors 
in PMI, as implied by Theorem |2.2| So, we only report the results for SN. 

For comparison, GloVe and two variants of word2vec (skip-gram and GBOW) vectors are trained. GloVe’s 
vectors are trained on the same co-occurrence as SN with the default parameter values]^ word2vec vectors 
are trained using a window size of 10, with other parameters set to default valuesj^ 


® http://nl p.sta nford .ed u / projects/glove / 
^ https: //code.google. com/p/word2vec/ 
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Natural logarithm of frequency 


Figure 2: The linear relationship between the squared norms of our word vectors and the logarithms of 
the word frequencies. Each dot in the plot corresponds to a word, where a;-axis is the natural logarithm of 
the word frequency, and y-axis is the squared norm of the word vector. The Pearson correlation coefficient 
between the two is 0.75, indicating a significant linear relationship, which strongly supports our mathematical 
prediction, that is, equation (2.4) of Theorem 


2.2 


5.1 Model verification 

Experiments were run to test our modeling assumptions. First, we tested two counter-intuitive properties: 
the concentration of the partition function for different discourse vectors c (see Theorem 2.1), and the 
random-like behavior of the matrix of word embeddings in terms of its singular values (see Theorem 4.1). 
For comparison we also tested these properties for word2vec and GloVe vectors, though they are trained by 
different objectives. Finally, we tested the linear relation between the squared norms of our word vectors 
and the logarithm of the word frequencies, as implied by Theorem |2.2| 


Partition function. Our theory predicts the counter-intuitive co ncen tration of the partition function 
Zf. — ^ random discourse vector c (see Lemma 2.1). This is verified empirically by 

picking a uniformly random direction, of norm ||c|| = 4//ru,, where is the average norm of the word 
vectors]^ Figure [^a) shows the histogram of Z^ for 1000 such randomly chosen c’s for our vectors. The 
values are concentrated, mostly in the range [0.9,1.1] times the mean. Concentration is also observed for 
other types of vectors, especially for GloVe and CBOW. 


Isotropy with respect to singular values. Our theoretical explanation of relations=lines assumes 
that the matrix of word vectors behaves like a random matrix with respect to the properties of singular 
values. In our embeddings, the quadratic mean of the singular values is 34.3, while the minimum non-zero 
singular value of our word vectors is 11. Therefore, the ratio between them is a small constant, consistent 
with our model. The ratios for GloVe, CBOW, and skip-gram are 1.4, 10.1, and 3.1, respectively, which are 
also small constants. 


Squared norms v.s. word frequencies. Figure shows a scatter plot for the squared norms of our 
vectors and the logarithms of the word frequencies. A linear relationship is observed (Pearson correlation 
0.75), thus supporting Theorem 2.2 The correlation is stronger for high frequency words, possibly because 


the corresponding terms have higher weights in the training objective. 


®Note that our model uses the inner products between the discourse vectors and word vectors, so it is invariant if the 
discourse vectors are scaled by s while the word vectors are scaled by 1/s for any s > 0. Therefore, one needs to choose the 
norm of c properly. We assume ||c||/ru, = '/dju ^ 4 for a constant k = 5 so that it gives a reasonable fit to the predicted 
dynamic range of word frequencies according to our theory; see model details in Section [2| 
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Relations 

SN 

GloVe 

CBOW 

skip-gram 


semantic 

0.84 

0.85 

0.79 

0.73 

G 

syntactic 

0.61 

0.65 

0.71 

0.68 


total 

0.71 

0.73 

0.74 

0.70 


adjective 

0.50 

0.56 

0.58 

0.58 


noun 

0.69 

0.70 

0.56 

0.58 


verb 

0.48 

0.53 

0.64 

0.56 


total 

0.53 

0.57 

0.62 

0.57 


Table 1: The accuracy on two word analogy task testbeds: G (the GOOGLE testbed); M (the MSR testbed). 
Performance is close to the state of the art despite using a generative model with provable properties. 

This correlation is much weaker for other types of word embeddings. This is possibly because they have 
more free parameters (“knobs to turn”), which imbue the embeddings with other properties. This can also 
cause the difference in the concentration of the partition function for the two methods. 


5.2 Performance on analogy tasks 


We compare the performance of our word vectors on analogy tasks, specifically the two testbeds GOOGLE 
and MSR (Mikolov et al., 2013a|c ). The former contains 7874 semantic questions such as “man:woman::king:7T\ 
and 10167 syntactic ones such as “run:runs::walk:!77^ The latter has 8000 syntactic questions for adjectives, 
nouns, and verbs. 

To solve these tasks, we use linear algebraic queries}^ That is, first normalize the vectors to unit norm 
and then solve “a:b::c:77” by 

argmin ||ua - vt - Vc + Vd\\l ■ 


(5.1) 


The algorithm succeeds if the best d happens to be correct. 

The performance of different methods is presented in Table Our vectors achieve performance compara¬ 
ble to the state of art on semantic analogies (similar accuracy as GloVe, better than word2vec). On syntactic 
tasks, they achieve accuracy 0.04 lower than GloVe and skip-gram, while GROW typically outperforms the 
others The reason is probably that our model ignores local word order, whereas the other models capture 
it to some extent. For example, a word “she” can affect the context by a lot and determine if the next word 
is “thinks” rather than “think”. Incorporating such linguistic features in the model is left for future work. 


5.3 Verifying relations=lines 


The theory in Section [^predicts the existence of a direction for a relation, whereas earlier [Levy and Goldberg 
(2014a) had questioned if this phenomenon is real. The experiment uses the analogy testbed, where each 
relation is tested using 20 or more analogies. For each relation, we take the set of vectors Vab = Va — Vb 
where the word pair (a, b) satisfies the relation. Then calculate the top singular vectors of the matrix formed 
by these fat’s, and compute the cosine similarity (i.e., normalized inner product) of individual Vab to the 
singular vectors. We observed that most (va — Vb)’s are correlated with the first singular vector, but have 
inner products around 0 with the second singular vector. Over all relations, the average projection on the 
first singular vector is 0.51 (semantic: 0.58; syntactic: 0.46), and the average on the second singular vector 
is 0.035. For example. Table shows the mean similarities and standard deviations on the first and second 
singular vectors for 4 relations. Similar results are also obtained for word embedings by GloVe and word2vec. 
Therefore, the first singular vector can be taken as the direction associated with this relation, while the other 
components are like random noise, in line with our model. 


®One can instead use the 3COSMUL in ( [Levy and Goldbergj [2014a| |, which increases the accuracy by about 3%. But it is 

not linear while our focus here is the linear algebraic structure. _ 

^®It was earlier reported that skip-gram outperforms CBOW i jMikolov et al.||2013a| [Pennington et aL||2014t . This may be 
due to the different training data sets and hyperparameters used. 
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relation 

1 

2 

3 

4 

5 

6 

7 

1st 

2nd 

0.65 ± 0.07 
0.02 ± 0.28 

0.61 ± 0.09 
0.00 ± 0.23 

0.52 ± 0.08 
0.05 ± 0.30 

0.54 ± 0.18 
0.06 ± 0.27 

0.60 ± 0.21 
0.01 ± 0.24 

0.35 ± 0.17 
0.07 ± 0.24 

0.42 ± 0.16 

0.01 ± 0.25 

relation 

8 

9 

10 

11 

12 

13 

14 

1st 

2nd 

0.56 ± 0.09 
0.00 ± 0.22 

0.53 ± 0.08 
0.01 ± 0.26 

0.37 ± 0.11 
0.02 ± 0.20 

0.72 ± 0.10 
0.01 ± 0.24 

0.37 ± 0.14 
0.07 ± 0.26 

0.40 ± 0.19 
0.07 ± 0.23 

0.43 ± 0.14 

0.09 ± 0.23 


Table 2: The verification of relation directions on 2 semantic and 2 syntactic relations in the GOOGLE 
testbed. Relations include cap-corn: capital-common-countries; cap-wor: capital-world; adj-adv: graml- 
adjective-to-adverb; opp: gram2-opposite. For each relation, take Vab = Va — Vb for pairs (a, b) in the 
relation, and then calculate the top singular vectors of the matrix formed by these VabS. The row with label 
“1st”/“2nd” shows the cosine similarities of individual Vab to the lst/2nd singular vector (the mean and 
standard deviation). 



SN 

GloVe 

CBOW 

skip-gram 

w/o RD 

0.71 

0.73 

0.74 

0.70 

RD(fc = 20) 

0.74 

0.77 

0.79 

0.75 

RD(fc = 30) 

0.79 

0.80 

0.82 

0.80 

RD(fc = 40) 

0.76 

0.80 

0.80 

0.77 


Table 3: The accuracy of the RD algorithm (i.e., the cheater method) on the GOOGLE testbed. The RD 
algorithm is described in the text. For comparison, the row “w/o RD” shows the accuracy of the old method 
without using RD. 


Cheating solver for analogy testbeds. The above linear structure suggests a better (but cheating) way 
to solve the analogy task. This uses the fact that the same semantic relationship (e.g., masculine-feminine, 
singular-plural) is tested many times in the testbed. If a relation R is represented by a direction fiji then the 
cheating algorithm can learn this direction (via rank 1 SVD) after seeing a few examples of the relationship. 
Then use the following method of solving “a:6::c:??”: look for a word d such that Vc — Vd has the largest 
projection on the relation direction for (a, 5). This can boost success rates by about 10%. 

The testbed can try to combat such cheating by giving analogy questions in a random order. But the 
cheating algorithm can just cluster the presented analogies to learn which of them are in the same relation. 
Thus the final algorithm, named analogy solver with relation direction (RD), is: take all vectors Va — Vb for 
all the word pairs (a, b) presented among the analogy questions and do /c-means clustering on them; for each 
{a,b), estimate the relation direction by taking the first singular vector of its cluster, and substitute that 
for Va — Vb in (5.11 when solving the analogy. Table shows the performance on GOOGLE with different 
values of k; e.g. using our SN vectors and A: = 30 leads to 0.79 accuracy. Thus future designers of analogy 
testbeds should remember not to test the same relationship too many times! This still leaves other ways to 
cheat, such as learning the directions for interesting semantic relations from other collections of analogies. 


Non-cheating solver for analogy testbeds. Now we show that even if a relationship is tested only 
once in the testbed, there is a way to use the above structure. Given “a:6::c:??,” the solver first finds the 
top 300 nearest neighbors of a and those of b, and then finds among these neighbors the top k pairs (o', b') 
so that the cosine similarities between Va' — Vb' and Va — Vb are largest. Finally, the solver uses these pairs 
to estimate the relation direction (via rank 1 SVD), and substitute this (corrected) estimate for Va — Vb in 
(5.1) when solving the analogy. This algorithm is named analogy solver with relation direction by nearest 
neighbors (RD-nn). Table shows its performance, which consistently improves over the old method by 
about 3%. 


15 






















SN 

GloVe 

CBOW 

skip-gram 

w/o RD-nn 

0.71 

0.73 

0.74 

0.70 

RD-nn {k = 10) 

0.71 

0.74 

0.77 

0.73 

RD-nn {k = 20) 

0.72 

0.75 

0.77 

0.74 

RD-nn (fc = 30) 

0.73 

0.76 

0.78 

0.74 


Table 4: The accuracy of the RD-nn algorithm on the GOOGLE testbed. The algorithm is described in the 
text. For comparison, the row “w/o RD-nn” shows the accuracy of the old method without using RD-nn. 


6 Conclusions 

A simple generative model has been introduced to explain the classical PMI based word embedding models, 
as well as recent variants involving energy-based models and matrix factorization. The model yields an 
optimization objective with essentially “no knobs to turn”, yet the embeddings lead to good performance 
on analogy tasks, and fit other predictions of our generative model. A model with fewer knobs to turn 
should be seen as a better scientific explanation {Occam’s razor), and certainly makes the embeddings more 
interpretable. 

The spatial isotropy of word vectors is both an assumption in our model, and also a new empirical 
finding of our paper. We feel it may help with further development of language models. It is important for 
explaining the success of solving analogies via low dimensional vectors (relations=LINES). It also implies 
that semantic relationships among words manifest themselves as special directions among word embeddings 
(Section]^, which lead to a cheater algorithm for solving analogy testbeds. 

Our model is tailored to capturing semantic similarity, more akin to a log-linear dynamic topic model. 
In particular, local word order is unimportant. Designing similar generative models (with provable and 
interpretable properties) with linguistic features is left for future work. 
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A Proofs of Theorem 1 


In this section we prove Theorem 2.2 and Lemma 2.1 (restated below). 

Theorem 


2.2 


Suppose the word vectors satisfy equation (|2.2[), and window size q = 2. Then, 


logp{w, w') = 
log p{w) = 


Iku. +Vw'\\i 

2d 

I l|2 


-2\ogZ±e, 


2d 


— log Z ±€. 


for e = 0{ez) + 0{l/d) + 0(e2)- Jointly these imply: 


PMI {w,w') = ± 0(e). 


(A.l) 

(A.2) 

(A.3) 


Lemma |2.1[ If the word vectors satisfy the bayesian prior v = s ■ v, where v is from the spherical Gaussian 
distribution, and s is a scalar random variable, then with high probability the entire ensemble of word vectors 
satisfies that 

Pr [(1 - e,)Z <Z,<{1 + ez)Z] >1-6, (A.4) 


c~C 


for Cz = 0{l/^/n), and 6 = exp(—n(log^ n)). 


We first prove Theorem |2.2| using Lemma |2.1[ and Lemma |2.1| will be proved in Section |A.1| Please see 
Section|^of the main paper for the intuition of the proof and a cleaner sketch without too many technicalities. 


Proof of Theorem 2.2 Let c be the hidden discourse that determines the probability of word w, and c' be 
the next one that determines w'. We use p(c'|c) to denote the Markov kernel (transition matrix) of the 
Markov chain. Let C be the stationary distribution of discourse vector c, and T be the joint distribution of 
(c, c'). We marginalize over the contexts c,c' and then use the independence of w,w' conditioned on c, c', 


p{w, w') = E 

(c,c')~'D 


exp((uu,,c)) exp((u„/,c')) 


(A.5) 


Zq Zc' 

We first get rid of the partition function Zc using Lemma [2T] As sketched in the main paper, essentially 


we will replace Z^ by Z in equation (A.5), though a very careful control of the approximation error is 


required. Formally, Let be the event that c satisfies 

(1 — ez)Z ^ Zc ^ {1 + ez)Z . 


(A.6) 


Similarly, let J ^2 be the even that c' satisfies {1 — €z)Z < Z^,' < (1 + e^).^, and let = J^ir\J^ 2 , and be its 
negation. Moreover, let Ijr be the indicator function for the event iF. Therefore by Lemma |2.1| and union 
bound, we have E[ljr] = Pr[J^] > 1 — exp(—r2(log^ n)). 

We first decompose the integral (A.5) into the two parts according to whether event F happens. 


p{w, w') = E 

(c,c')~D 

+ E 


1 


■ exp((uu,, c)) exp((uu,/, c'))ljF 


{c,c')r^'D _ZcZ( 


■ exp ((uu,, c)) exp( (uu,', c')) 1 


We bound the first quantity on the right hand side using (2.2) and the definition of F. 

1 


E 

(c,c')~X) 


exp((z;„, c)) exp((uu,/, c'))l 




1 


< (l + Czf^ E _ [exp((u^,,c))exp((uu,/,c'))l^] 


(c.c')~D 


(A.7) 


(A.8) 
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For the second quantity of the right hand side of (A.7|, we have by Cauchy-Schwartz, 


E 

(c,c')~X) 


ZnZ, 


< 


E 

(c,c')~D 


1 


< E 


1 

Using the fact that Zc > 1, then we have that 


exp ((uu,, c)) exp( (uu,/ , c )) 

exp((uu,,c))^ly^ 

DC 

E 

c')~D 

E [1^] 

c'\c 

)(? 



Jc,d ^d 


J \\‘2 


E 


^exp((z;^,c))2 E [ly^] 

c'|c 


< E 


c\c' 


exp((wu,,c))^ E [ly^] 


We can split that expectation as 


E 


exp((riu,,c))^l(^^^c)>o E [1 


'k 




■E 


^Xp(('C-u;, c)) l(t;^,c )<0 E [Ijr] 

c'jc 


The second term of (A. 10) is upper bounded by 


E [Ijr] < exp(-U(log^n)) 


(A.9) 


(A.IO) 


We proceed to the first term of (|A.10 ) and observe the following property of it: 


E 


exp((z)u,,c))^l(„„,c)>o E [ly^] 

c'\c 


<E 


exp((a?;u,,c))^l(„„,c)>o E [ly^] 

d\c 


< E 


exp((«■(;„,, c))^ E [ly^] 

c'|c 


where a > 1. Therefore, it’s sufficient to bound 


E 


exp((uu,,c))^ E [Ijr] 

c'\c 


when ||uu,|| = fl{Vd). 

Let’s denote by z the random variable 2 {vw,c). 

Let’s denote r(z) = Ec'|z[ly^] which is a function of z between [0,1]. We wish to upper bound Ec [exp(z)r(z)]. 
The worst-case r(z) can be quantified using a continuous version of Abel’s inequality as proven in Lemma [A.4[ 
which gives 

E[exp(z)r(z)] < E [exp(z)l[t^ +oo](^)] (A.ll) 

C 

where t satisfies that Ec[l[t,+oo]] = Pr[-z ^ = Ec[?'(2:)] < exp(—r2(log^ n)). Then, we claim Pr[z > t] < 

exp(—U(log^ n)) implies that t > U(log'® n). 

If c were distributed as A/'(0, ^F), this would be a simple tail bound. However, as c is distributed uniformly 
on the sphere, this requires special care, and the claim follows by applying Lemma instead. Finally, 
applying Corollary |A.3[ we have: 

E[exp(z)r(z)] < E [exp(z)l[t,+oo](z)] = exp(-U(log^'® n)) (A.12) 

We have the same bound for c' as well. Hence, for the second quantity of the right hand side of ( |A7t , 
we have 


E 

(c,c')~D 


■ exp((u„, c)) exp((ui„/, c'))ly' 


1 


— exp((uu„c)) E [ly^] 

c'|c 


< exp(—H(log^ ® n)) 


1/2 


E 


^exp((u^/,c')) E [ly^] 




c\d 


1/2 


(A.13) 
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where the first inequality follows from Cauchy-Schwartz, and the second from the calculation above. Com¬ 


bining (A.71, (A.8| and (A.131, we obtain 


p{w,w') < (1-f £ 2 )^^ E [exp((r;u,,c))exp((r;^/,c'))lj.] -f ^ exp(-fl(log^ ® n)) 

^ (c,c')~X> Tl 

<(l + ez)^^, E _ [exp((r;u,,c))exp(('!;u,/,c'))]-f (5o 


(c,c')~V 


where So = exp(—f2(log^ ® n))Z^ < exp(—fl(log^'® n)) by the fact that Z < exp(2K)n = 0{n). Note that 
K is treated as an absolute constant throughout the paper. On the other hand, we can lowerbound similarly 

p{w,w') > {1 - E [exp(('Uu,,c))exp((u„/,c'))ly^] 

^ {c,C')r^V 

>(l-ez)^^ E [exp((i;u,,c))exp((u„/,c'))] - ^exp(-fl(log^'®n)) 

(c,c')~-D 

>(l-ez)^^ E [exp((i;u,, c))exp( (?;„/, c'))] - (5o 

Z‘‘ {c,c')~V 

Taking logarithm, the multiplicative error translates to a additive error 

logp(u;, w') = log ( E [exp((i;u,,c))exp((?;u,',c'))] ± ^0 ) - 2 logZ-f 2log(l ± £ 2 ) 


For the purpose of exploiting the fact that c, c' should be close to each other, we further rewrite logp(w, w') 
by re-organizing the expectations above. 


logp(ui,'u;') = log ( E 


± (5o - 2 log Z -f 2 log(l ± £^) 


exp((z;^,,c)) E [exp((r)u,/,c'))] 


= log E [exp((z)u,, c))A(c)] ± (5o I - 2 log Z + 2 log(l ± e^) 


(A.14) 


where the inner integral which is denoted by A(c), 


A(c) := E [exp((z;^,/,c'))] 

c'|c 


Since ||f,u|| < K'/d. Therefore we have that {vyj,c — c') < ||r;uj||||c — c'|| < K'</d\\c — c'\\. 
Then we can bound A(c) by 


A(c) = E [exp((z;u,',c'))] 

c'|c 

= exp((ri„/,c)) E [exp((uu,/, c'- c))] 

c'|c 

< exp((r!u,/, c)) E [exp(Kv^||c - c'||)] 

c'|c 

< (l-l-£2)exp((nu,',c)) 

where the last inequality follows from our model assumptions. To derive a lower bound of A(c), observe that 
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E [exp(K-\/d||c — c'll)] + E [exp(—/c-\/d||c — c'||)] > 2 

c'|c c'|c 


Therefore, our model assumptions imply that 


Hence, 


E [exp(-KV^||c - c'll)] > 1 - e 2 

c'|c 


H(c) 


Therefore, we obtain that A{c) = 
the equation (A. 14), we get that 


= exp((uu,/, c)) E exp((uu,/, c'- c)) 

c'\c 

> exp((uu,', c)) E exp(-KV^||c - c'll) 

c'\c 

> (1 - e2)exp((c„',c)) 

(1 ± £ 2 ) exp((uuj', c)). Plugging the just obtained estimate of A{c) into 


logp('u;,ui') = log I^E [exp((u^,,c))A(c)] ± Joj - 21ogZ + 21og(l ± e^) 

= log ^E [(1 ± £ 2 ) exp((u„, c)) exp((z)u,/, c))] ± i5o^ - 2 log Z + 2 log(l ± £^) 
= log (e [exp((u^ +v^>,c))] ± ho ) - 2logZ + 2log(l ± £^) + log(l ± £ 2 ) 


(A.15) 


Now it suffices to compute Ec[exp((?;u, + z;u,/,c))]. Note that if c had the distribution Af{0, g/), which 
is very similar to uniform distribution over the sphere, then we could get straightforwardly Ec[exp((uu, + 
^^lu'jc))] = exp(||cu, + Viiji /{2d)). For c having a uniform distribution over the sphere, by Lemma A.5[ the 
same equality holds approximately, 

E[exp((ui„+ u„/,c))] = (l±£ 3 )exp(||ui„ + ?;„/|p/( 2 d)) (A.16) 

C 

where £3 = 0{l/d). 

Plugging in equation ( |A.16 ) into equation (A.15), we have that 

logp('u;, w') = log ((1 ± £ 3 )exp(||?;„ + Uuj'|p/(2d)) ± ho) - 2 log Z + 2 log(l ± e^) + log(l ± £ 2 ) 

= ll'i'tu + Vw' |P/(2d) + 0 (^ 3 ) + 0{6q) — 2 log Z ± 2ez ± £2 


where hg = ho • (Ec~c[exp((?;u, + u^,/, c))]) ^ = exp(-H(log^ ® n)). Note that £3 = 0{l/d), £z = 0{l/y/n), 
and £2 by assumption, therefore we obtain that 

logp(w,'u;') = Ti'yllUu, Tu^'lP - 21 ogZ± 0 (£^) + 0 (£ 2 ) + d{l/d). 


□ 

The following lemmas are helper lemmas that were used in the proof above. We use to denote the 
uniform distribution over the unit sphere in 

Lemma A.l (Tail bound for spherical distribution). If c ^ C‘^, v £ is a vector with ||w|| = Vt{^fd) and 
t = w(l), the random variable z = {v,c) satisfies Pt[z > t ] = \ 
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Proof. If c = (ci,C 2 ,... ,Cd) ~ C‘^, c is in distribution equal to ■ ■ •) where the Ci are i.i.d. 

samples from a univariate Gaussian with mean 0 and variance By spherical symmetry, we may assume 
that V = (||?;||, 0,..., 0). Let’s introduce the random variable r = Yl'i =2 ■ Since 


Pr [{v, c) > t] = Pr 




< Pr 


IzjIIci . 1 

5 - ' - 2 


Pr 


r > 


Pr 


k||ci I 1 

II5II - ' - 2 


Pr 


r > 


it’s sufficient to lower bound Pr [r < 100] and Pr ||'y||]^ t \ r < 100 . The former probability is easily 
seen to be lower bounded by a constant by a Chernoff bound. Consider the latter one next. It holds that 


Pr 


kllpn- > t I r < 100 


= Pr 

IV 




f^ ■ r 




r < 100 


> Pr 


Cl > 


100^2 

U|| 2 -t 2 


Denoting i = , by a well-known Gaussian tail bound it follows that 




where the last equality holds since ||u|| = U{'/d) and t = w(l). □ 

Lemma A.2. If , v G is a vector with ||n|| = 0(-\/d) and t = w(l), the random variable z = {v,c) 

satisfies 

E [exp( 0 )l[t_+oo]( 2 :)] = exp(-fl(t^)) -h exp(-n(d)) 


Proof. Similarly as in Lemma 


A.l 


if c = (ci,C 2 ,... ,Cd) ~ C"*, c is in distribution equal to (^|||||, |j||, ■ ■ ■, p|) 
where the Ci are i.i.d. samples from a univariate Gaussian with mean 0 and variance g. Again, by spherical 
symmetry, we may assume v = (||u||, 0,..., 0). Let’s introduce the random variable r = Then, for 

an arbitrary u > 1, some algebraic manipulation shows 

Pr [exp((i;,c)) l[ 4 ,+oo]((?^, c)) > u] = Pr [exp ((u, c)) >uA{v,c)>i\ = 


Pr 


exp \\v 


Cl 


> uA Iklljrli > u 


= Pr 


Cl = max 



(A.17) 


where we denote u = logu. Since ci is a mean 0 univariate Gaussian with variance and ||i;|| = fl^Vd) we 
have Va: G K 


Pr 


Cl > 


= 0 e 




Next, we show that r is lower bounded by a constant with probability 1 — exp(—n(d)). Indeed, r is in 
distribution equal to ^Xd-i^ where Xk i® ^ Chi-squared distribution with k degrees of freedom. Standard 

concentration bounds (?) imply that > 0,Pr[r — 1 < < exp(—^). Taking f = adior a & constant 

implies that with probability 1 — exp(—n(d)), r > M for some constant M. We can now rewrite 


Pr 


Cl > 
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Pr 


Cl > 


- r > M 

\\vP-x^ ' - 


Pr[r > M] + Pr 


Cl > 


- r < M 

||i;|P-x 2 I - 


Pr[r < M] 


The first term is clearly bounded by e and the second by exp(—n((i)). Therefore, 


Pr 


Cl > 


= O (max (exp (—17 (x^)) ,exp (—17 (d)))) 


Putting [A. 17| and | A. 18| together, we get that 

Pr [exp {{v, c)) l[t_+oo]((^^, c)) > u] = O ^max ^exp ^—17 ^min (^d, (max (it, 7))' 
(where again, we denote u = log u) 

For any random variable X which has non-negative support, it’s easy to check that 


(A.18) 


(A.19) 


poo 

E[A] = / Pr[A > x]d 2 

Jo 


Hence, 


roo rexp(||«||) 

E [exp(z)l[t^+oo](^)] = / Pr [exp( 0 )l[t,+oo] (z) >u]du= I Pr [exp(z)l[t_+oo]( 2 :) > u] du 

Jo Jo 

To bound this integral, we split into the following two cases: 

• Case > d: max (12,7) > 7, so min ^d, (max (12,7))^^ = d. Hence, 


A.19 


implies 


E [exp(z)l[t,+oo](2:)] = exp(||i;||)exp(-17(d)) = exp(-17(d)) 
where the last inequality follows since ||i;|| = 0{'/d). 

Case 7^ < d: In the second case, we will split the integral into two portions: u G [0, exp(7)] and 
u G [exp(7),exp(||i;||)]. 

When u G [0,exp(7)], max (12, 7) = 7, so min(d, (max (12, 7))^) = 7^. Hence, 

/■exp(t) 

/ Pr [exp( 2 ;)l[( _|_oo](z) > u] dit < exp(7) exp(—17(7^)) = — exp(17(7^)) 

Jo 

When u G [exp(7), exp(||i;||)], max (12, 7) = u. But u < log(exp(||i;||)) = 0{Vd), so min(d, (max (12, 7))^) = 
u. Hence, 

rexp(lhll) cexp(||t;||) 

/ Pr [exp(z)l[t_+oo](- 2 ) > u] dw < / exp(-(log(u))^)dM 

^ exD(t) ^ exx}(t) 


Making the change of variable u = log(ii), the we can rewrite the last integral as 

rlhll 

/ exp(— 12 ^) exp{u)du = 0 (exp(—7^)) 


where the last inequality is the usual Gaussian tail bound. 
In either case, we get that 

^exp(||u|r 


ne^Pdlvw) 

/ Pr [exp(z)l[t_+oo]( 2 :) >u]du = exp(-I7(7^)) -|- exp(-I7(d))) 

Jo 


which is what we want. 


□ 
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As a corollary to the above lemma, we get the following: 

Corollary A.3. If c ^ , v is a vector with ||z;|| = Q{\/d) and t = n(log'® n) then 

E [exp(z)l[t_+^](z)] = exp(-fl(log^-®n)) 


Proof. We claim the proof is trivial if d = oflog^ n). Indeed, in this case, exp((ti, c)) < exp(||n||) = 
exp{0{Vd)). Hence, 


E [exp(z)l[t^+^](z)] = exp(0(yd))E[l[t,+oo](2:)] = exp(0(yd)) Pr[z > t] 


Since by Lemma A.l[ Pr[z > t] < exp(—H(log n), we get 


E [exp(2;)l[t,+oo](z)] = exp(0(v^) - H(log^n)) = exp(-H(log^ ® n)) 


as we wanted. 

So, we may, without loss of generality assume that d = f2(log^ n). In this case, Lemma 


A.2 


implies 


E [exp{z)l[t_+^] (x)] = exp(- log^ ® n) + exp(-H(d))) = exp(- log^'® n) 
where the last inequality holds because d = H(log^ n) and t^ = H(log'® n), so we get the claim we wanted. 


□ 


Lemma A.4 (Continuous Abel’s Inequality). Let 0 < r{x) < 1 be a function such that such that E[r(a:)] = p. 
Moreover, suppose increasing function u{x) satisfies that E[|u(a:)|] < oo. Let t be the real number such that 
E[l[t,+oo]] = P- Then we have 

E[M(a;)r(x)] < E[u(a:)l[t_+oo]] (A.20) 

Proof. Let G{z) = f{x)r{x)dx, and H{z) = f{x)l[t^.^.ao]ix)dx. Then we have that G{z) < H{z) for 
all z. Indeed, for z >t, this is trivial since r{z) < 1. For z < t, we have H{z) = E[l[t,+oo]] = P = E[r(a:)] > 
f{x)r{x)dx. Then by integration by parts we have, 


/ OO poo 

u{x) f {x)r{x)dx = — u{x)dG 

-OO J —oo 

/ + 00 

G{x)u' {x)da 

-OO 
XD 

H{x)u' {x)dx 

) 

u[x)f{x)l^t,+oc,]{x)dx, 


P + OO 


where at the third line we use the fact that u{x)G{x) —)■ 0 as a; —>■ oo and that u'{x) > 0, and at the last line 
we integrate by parts again. □ 

Lemma A.5. Let v G be a fixed vector with norm ||u|| < k\/ d for absolute constant k. Then for random 
variable c with uniform distribution over the sphere, we have that 

logE[exp((u, d))] = ||u||^/2d±e (A.21) 


where e = 0 ( 3 ). 
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Proof. Let g S J\f{0,I), then g/\\g\\ has the same distribution as c. Let r = ||u||. Since c is spherically 
symmetric, we could, we can assume without loss of generality that v = (r, 0, Let x = gi and 

y = V9i + --- + g^. Therefore x € Af(0, 1) and has distribution with mean d — 1 and variance 0{d). 

Let be the event that a; < 20 log d and 1.5-\/d > y > O.bVd. Note that the Pr[J^] > 1—exp(—f2(log^ ®(d)). 
By Proposition]^ we have that E[exp((u, c))] = E[exp((u, c)) | • (1 ± n{— log^'® d)). 

Conditioned on event J-, we have 


E[exp((u, c)) I J"] = E 


exp( 


rx 

y/x"^ + y'^ 


) 1-^ 


= E 

= E 

= E 


,rx 
exp(- 


y y^Jxf+^{y + ^/xF+^) 
rx 


) 1-^ 


exp( —) • exp(— - ) I T 

y y^Jx^ + y^[y P ^Jx^ + y^) 


exp( —) I F 
V 


(l±0(f^)) 


(A.22) 


where we used the fact that r < K\fd. Let £ be the event that l.bVd > y > 0.5Vd. By using Proposition]^ 
we have that 

E[exp{rx/y) \ F] = E[exp(rx/y) \ £] ± exp(-f2(log^(d)) (A.23) 

Then let z = y'^j{d— 1) and w = z — 1. Therefore z has distribution with mean 1 and variance l/(d— 1), 
and w has mean 0 and variance l/(d — 1). 


E 


exp( —) I £ 

y 


= E[E[exp(ra;/j/) \y\\£]= E[exp(r^/?/2) | 

= E[exp(rV(d- 1) • l/z^) \ £] 

Oqn -I— 7/;^ 

= E[exp(rV(d-l)-(l + ^Y^))l^] 

O7/1 _|_ 7/1^ 

= exp(rV(d - 1)) E[exp(l + I £] 

= exp(r^/ {d — 1)) E[1 + 2ri; ± 0{w^) \ £] 

= exp(rV(d-l)^)(l±l/d) 


where the second-to-last line uses the fact that conditioned on 1/2 > £, ru > —1/2 and therefore the Taylor 
expansion approximates the exponential accurately, and the last line uses the fact that | E[rc I '^^]| = 0{l/d) 
and E[w^ | £] < 0{l/d). Combining the series of approximations above completes the proof. 


□ 


We finally provide the proofs for a few helper propositions on conditional probabilities for high probability 
events used in the lemma above. 

Proposition 1. Suppose x ^A/’(0,cr^) witha = 0(1). Then for any event£ withPT[£] = 1 —0(—n(log^ d)), 
we have that E[exp(a;)] = E[exp(a;) | £] ± exp(—f2(log^(d)). 

Proof Let’s denote by £ the complement of the event £. We will consider the upper and lower bound 
separately. Since 

E[exp(x)] = E[exp(a;) | £] Pr[£l] + E[exp(x) | £] Pr[f] 

we have that 

E[exp(a;)] < E[exp(a;) | £] + E[exp(x) | £] Pr[f] (A.24) 
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and 


E[exp(a;)] > E[exp(a;) | f](l — exp(— r2(log^ d))) > E[exp(x) | £] — E[exp(a;) | £] exp(— r2(log^ d))) (A.25) 


Consider the upper bound (A.241 first. To show the statement of the lemma, it suffices to bound 
E[exp(a;) | £]Vv[£]. 

Working towards that, notice that 

E[exp(a;) | £] Pr[£] = E[exp(a:)l£] = E[exp(x) E[l£|a;]] = E[exp(x)r(x)] 
if we denote r[x) = E[l£|a;]. We wish to upper bound E[exp(x)r(a;)]. By Lemma [W4 we have 

E[exp(x)r(x)] < E[exp(a:)l[t_oo]] 

where t is such that E[l[t,oo]] = E[r'(a:)]. However, since E[r'(a;)] = Pr[£] = exp(—f2(log^ d)), it must be the 
case that t = H(logd) by the standard Gaussian tail bound, and the assumption that a = 0(1). In turn, 
this means 


E[exp(a;)l[t^oo]] < 


1 




e e 


T\p2/K . 


e (cr -i) i dx = e 


2 1 


^+oo 




e 2^^ dec' 


't/o 


where the last equality follows from the change of variables x = ax'. However, 


^ + O0 




e ^^^dx' 


tjo 


is nothing more than Pr[a;' > ^], where x' is distributed like a univariate gaussian with mean ^ and variance 
1. Bearing in mind that cr = 0(1) 


.2 1 


e * 




/ da;'= exp(—f2(t^)) = exp(—0(log^ d)) 

J tj<7 


by the usual Gaussian tail bounds, which proves the lower bound we need. 

We proceed to consider the lower bound |A.25[ To show the statement of the lemma, we will bound 
E[exp(a;) | £]. Notice trivially that since exp(a;) > 0, 


E[exp(a;) | £] < 


E[exp(a;)] 


Since Pr[f] > 1 — exp(H(log^ d)), < 1 + exp(0(log^)). So, it suffices to bound E[exp(a;)]. However, 


E[exp(x)] = 


f* + 00 


G's/ZTT Jt— — c 


ri/27r Jt=- 


+ 00 , 1 

e <T 2 ) + 4 . dx = 

00 


^ + 00 


Jt=-c 


3 (P 2)^+"4 dx' 


where the last equality follows from the same change of variables x = ax' as before. Since e 2 )^da;' = 

we get 


^-j-OO 2 2 

/ = gT = 0(1) 

J t— — oo 


A.25 

-□ 


Jt=- 

Putting together with the estimate of we get that E[exp(a:) | £] = 0(1). Plugging this back in 
we get the desired upper bound. 

Proposition 2. Suppose c ^ C and v is an arbitrary vector with ||u|| = 0{Vd). Then for any event £ with 
Pr[£] > 1 — exp(—H(log^ d)), we have t/iat E[exp((i;, c))] = E[exp((u, c)) | £] ± exp(—log^'® d). 
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Proof of Proposition^ Let z = {v, c). We proceed similarly as in the proof of Proposition]^ We have 

E[exp(z)] = E[exp(z) | £] Pr[£:] + E[exp(z) | £] Pr[£] 

and 

E[exp(z)] < E[exp(z) | £] + E[exp(z) | £] Pr[£] (A.26) 

and 

E[exp(z)] > E[exp(z) | £] Pr[£^] = E[exp(z) | £] — E[exp(z) | £] exp(—n(log^ d)) (A.27) 

We again proceed by separating the upper and lower bound. 

We first consider the upper bound |AT^ 

Notice that that 

E[exp(z) I f]Pr[5] =E[exp(z)l£] 

We can split the last expression as 

E [exp((u^,c))l(j,^,c)>o%] +E [exp((uu,,c))l^„„,c)<o%] ■ 

The second term is upper bounded by 

E)!^] < exp(—(log^ n)) 


We proceed to the first term of (A.IO) and observe the following property of it: 

E [exp(c))l(^^,c)>o If] < E [exp((afi„, c))l(„^,c)>o If] < E [exp((aUu,, c))ly] 
where a > 1. Therefore, it’s sufficient to bound 

E[exp(z)%] 

when lluu,!! = Q{Vd). Let’s denote r(z) = E[lf |z:]. 

Using Lemma [A.4[ we have that 


(A.28) 


E[exp(z)r(z)] < E [exp(z)l[t_+^](z)] 

C 

where t satisfies that Ec[l[t,+oo]] = Pr[-z > = Ec[»"(-z)] < exp(—r2(log^ d)). Then, we claim Pr[z > t] < 

exp(—n(log^ d)) implies that t > U(log ® d). 

Indeed, this follows by directly applying Lemma o Afterward, applying Lemma |A.2[ we have: 

E[exp(z)r(z)] < E [exp(z)l[t_+oo](^)] = exp(-n(log^-® d)) (A.29) 

which proves the upper bound we want. 

We now proceed to the lower bound |A.27| which is again similar to the lower bound in the proof of 
Proposition]^ we just need to bound E[exp(z) | £]. Same as in Proposition]^ since exp(a;) > 0, 

E[exp(z)] 


E[exp(z) I £] < 


Pr[f] 


Consider the event £' : z < t, for t = 0(log ® d), which by Lemma A.l satisfies Pr[f'] > 1 — exp(—r2(log^ d)). 
By the upper bound we just showed, 

E[exp(z)] < E[exp(z) | £'] + exp(-fI(log^ n)) = 0(exp(log ® d)) 

where the last equality follows since conditioned on z = 0(log'® d). Finally, this implies 


E[exp(z) I £] < ^j^^O(exp(log'®d)) = 0(exp(log ® d)) 


where the last equality follows since Pr[f] > 1 — exp(—n(log^ n)). Putting this together with A.27 
E[exp(z)] > E[exp(z) | £] Pr[f ] = E[exp(z) | £] — E[exp(z) | £\ exp(—n(log^ d)) > 
E[exp(z) I £] — 0(exp(log'® d)) exp(—U(log^ d)) > E[exp(z) | £] — exp(—U(log^ d)) 
which is what we needed. 


we get 


□ 
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A.l Analyzing partition function 

In this section, we prove LemmaWe basically first prove that for the means of are all (1 + o(l))-close 
to each other, and then prove that Z^ is concentrated around its mean. It turns out the concentration part 
is non trivial because the random variable of concern, exp((i;u,, c)) is not well-behaved in terms of the tail. 
Note that exp((uu,,c)) is NOT sub-gaussian for any variance proxy. This essentially disallows us to use an 
existing concentration inequality directly. We get around this issue by considering the truncated version of 
exp((uu,, c)), which is bounded, and have similar tail properties as the original one, in the regime that we 
are concerning. 

We bound the mean and variance of Z^^ Hrst in the Lemma below. 

Lemma A. 6 . For any fixed unit vector c S we have that E[^c] ^ ^ V[Zc] < Oin). 

Proof of Lemma \A.6[ Recall that by definition 

Zc = y^exp(('c^,,c)). 

W 

We fix context c and view v^^s as random variables throughout this proof. Recall that Vu, is composed of 
Vw = s-u! ■ v-u!, where s^j is the scaling and Vw is from spherical Gaussian with identity covariance Idxd- Let 
s be a random variable that has the same distribution as Sw 
We lowerbound the mean of Z^ as follows: 


E[^c] = «E [exp((uu,, c))] > nE [1 -I- {vw, c)] = n 

where the last equality holds because of the symmetry of the spherical Gaussian distibution. On the other 
hand, to upperbound the mean of Zc, we condition on the scaling s^,. 


E[Zc] = nE[exp((uu,,c))] 

= nE[E[exp((uu,,c)) | s^]] 

Note that conditioned on s^,, we have that {vw,c) is a Gaussian random variable with variance 
Therefore, 


E[exp((uu,,c)) I Su,] = 


^exp(-^)exp(*)<ix 

{x — (T^)^ 2 /r,N 

-I- (7^I2)dx 


= [ ^:^exp(- 
Jx crv27r 

= exp(cr^/ 2 ) 


2u2 


It follows that 

E[^c] = nE[exp(cr 2 / 2 )] = nE[exp(s^/ 2 )] = nE[exp(s^/ 2 )]. 
We calculate the variance of Zc as follows: 


V[Zc] = ^V[exp((uu,,c))] < n E[exp(2 (?;„,, c))] 

W 

= nE [E [exp(2(z)u,,c)) | s^,]] 

By a very similar calculation as above, using the fact that 2{vw,c) is a Gaussian random variable with 
variance dcr^ = 45 ^^ 


E [exp(2(w^,,c)) I Snj] = exp(2(T^) 
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Therefore, we have that 


V[Zc] < nE[E[exp(2(u„,c)) | s^,]] 

= nE [exp(2(T^)] = nE [exp(2s^)] < An 


for A = exp(8n^) a constant, and at the last step we used the facts that s < k a.s. 


□ 


Now we are ready to prove Lemma 2.1 


Proof of Lemma \2.1\ We fix the choice of c, and the proving the concentration using the randomness of Vw ’s 
first. Note that that exp((uu,,c)) is neither sub-Gaussian nor sub-exponential (actually the Orlicz norm of 
random variable exp((uu,,c)) is never bounded). This prevents us from applying the usual concentration 
inequalities. The proof deals with this issue in a slightly more specialized manner. 

Let’s define J'u, be the event that |(uu,,c)| < |logn. We claim that Pr[J'u,] > 1 —exp(—n(log^ n)). Indeed 
note that (fu,,c) | s™ has a Gaussian distribution with standard deviation Su)||c|| = < 2k a.s. Therefore 

by the Gaussianity of {vuj , c) we have that 


Pr[|(w^,,c)| > 7?logn I Su,] < 2exp(-fl(-log^ = exp(-r2(log^ n)). 


where r2(-) hides the dependency on k which is treated as absolute constants. Taking expectations over s^, 
we obtain that 

Pr[7'u,] = Pr[|(z;u,,c)| < ^logn] > 1 - exp(-fI(log^ n)). 

Note that by definition, we in particular have that conditioned on it holds that exp ((tin,, c)) < 

Let the random variable have the same distribution as exp((w„, c))|jr^. We prove that the random 
variable concentrates well. By convexity of the exponential function, we have that the mean 

of Z' is lowerbounded 


E[Z'] = uE[exp((uu,,c))|jr^] > nexp(E[(z;u,,c)|jr„]) =n 
and the variance is upperbounded by 


V[Z'] < nE [exp((n^,c))^|jF„] 

< j E [exp((u^,,c})^] 

< —p^;^An < l.lAn 

" Pr[J^^] - 

where the second line uses the fact that 


E [exp((?;i„,c))^] = Pr[J'u,] E [exp((uu,, c))^|JPu,] -k Pr[J^u,] E [exp(('(;„, c))^|J^uj] 

> Pr[7'u,]E [exp((nu,,c))^|J'u,] . 

Moreover, by definition, for any rc, 1X^,1 < y/n. Therefore by Bernstein’s inequality, we have that 


Pr [|Z' — E[^c]| > en] < exp(— 


le^n^ 


l.lAn -I- \^n ■ en 
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Note that E[Z'^] > n, therefore for e we have, 

Pr IIZ' - E[^']| > eE[^']] < Pr [jZ' - E[^']| > en] < exp(- % -) 

An + I 

< exp(—r2(min{e^n/A, ey^})) 

< exp(—0(log^ n)) 

Let = UyjJ^w be the union of all J^yj- Then by union bound, it holds that Pr[^] < < 

n ■ exp(—Sl(log^ n)) = exp(—n(log^ n)). We have that by definition, Z' has the same distribution as Zc|jr. 
Therefore, we have that 

Pr[|Z, - E[Z']| > eE[Z;] | < exp(-L!(log2 n)) (A.30) 

and therefore 


Pr[|Z, - E[Z']| > eE[Z']] = Pr[^] • Pr[|Ze - E[Z']| > eE^] | + Pr[^] Pr[|Z, - E[Z']| > eE[Z'] | JF] 

<Pr[|Z,-E[Z']| >eE[Z'] I J-]+Pr[jr] 

< exp(—f](log^ n)) (A.31) 


where at the last line we used the fact that Pr[j'] < exp(—n(log^ n)) and equation (A.30). 

Let Z = E[Z'] = E[exp((z)u,, c)) | | {vy,,c) \ < ilogn] (note that E[Z'] only depends on the norm of ||c|| 
which is equal to 1). Therefore we obtain that with high probability over the randomness of u^j’s, 

(l-ez)^<^c<(l + e.)^ (A.32) 

Taking expectation over the randomness of c, we have that 

Pr [ (A.32) holds] > 1 — exp(—0(log^ n)) 


Therefore by a standard averaging argument (using Markov inequality), we have 


Pr 

Vw L C 


Pr[ (A.321 holds] > 1 — exp(—n(log n)) > 1 — exp(—f2(log n)) 


For now on we fix a choice of so that Prc[ (A.32) holds] > 1 — exp(—r2(log^ n)) is true. Therefore in the 
rest of the proof, only c is random variable, and with probability 1 — exp(—n(log^ n)) over the randomness 
of c, it holds that, 

(1 - e,)Z < Z, < (1 + e,)Z. (A.33) 

□ 


B Maximum likelihood estimator for co-occurrence 

Let L be the corpus size, and Xy,^yji the number of times words w, w' co-occur within a context of size 10 
in the corpus. According to the model, the probability of this event at any particular time is logp(i/;,w') oc 
+ ftu'|l 2 • Successive samples from a random walk are not independent of course, but if the random walk 
mixes fairly quickly (and the mixing time of our random walk is related to the logarithm of the number 
of words) then the set of Xy,^y,Fs over all word pairs is distributed up to a very close approximation as a 
multinomial distribution Mul(L, {^(w, rc')}) where L = total number of word pairs in 

consideration (roughly lOL). 
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Assuming this approximation, we show below that the maximum likelihood values for the word vectors 
correspond to the following optimization, 


mm 

M,c 




(Objective SN) 


Now we give the derivation of the objective. According to the multinomial distribution, maximizing the 
likelihood of is equivalent to maximizing 


^ = log| j = ^ X^^^>logp{w,w'). 

\{w,w') j {w^w') 

To reason about the likelihood, denote the logarithm of the ratio between the expected count and the 
empirical count as 


X'ln^w' — log 


Lp{w, w') 


Note that 


^ X^^^i\ogp{w,w') 

{w,w') 


^ ^ ^W,'U 

{w,w') 


log- 


a:„ 


log 


Lp{w, w') 


X, 


w,w' 


x„ 


_ ^ 1 Xyj^y^' (Lp(w,w') 

— / ^ ^w^w' lC)§ ^ ^w,w' lC)§ 

{w^w') {w^w') 

— C “h ^ ^ ^w^w' ^w^w' 

{w^w') 


X 

where we let c denote the constant X!(u; w^) ^w,w' log ■ Furthermore, we have 

{w,w') 

= X X 

{w,w') 

= Alu,^u,/(1 + Au,,u,'+ + 0(|Au,,ju'|^)) 

(it;,It;') 

and also L = X)(u,,u;') Xw,w'- So 


(B.l) 


Xw,w'Xw,w'— — I A^^^,//2 + Xuj,u)'0(|Auj,u,' 1^) 

(it;,It;') 


v{w,w') 


(it;,It;') 


Plugging this into (3.2) leads to 


' — i— A^ ,„//2 + Auj,u,'0(|Auj^^i 


u,w') 


(it;,It;') 


(B.2) 
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When the last term is much smaller than the first term on the right hand side, maximizing the likelihood 
is approximately equivalent to minimizing the first term on the right hand side, which is our objective: 

[\\vw+Vwi\\l/{2d)+\ogL-2\ogZ^ 

{w^w') {w^w') 


where Z is the partition function. 


We now argue that the last term is much smaller than the first term on the right hand side in (B.2). For 


a large the is close to 0 and thus the induced approximation error is small. Small Xy^^yji’s only 

contribute a small fraction of the final objective (3.3), so we can safely ignore the errors. To see this, note 
that the objective Y.(w w')^w,w'^1, w' and the error term 1^) differ by a factor of 


|A^^^z| for each X^ ^n 


{w,w') ^W,w' 

'. For large ,„/’s, |Au,,u,'| 1, and thus their corresponding errors are much smaller 

than the objective. So we only need to conside r the Xyj yjds that are small constants. The co-occurrence 
counts obey a power law distribution (see, e.g. ( {Pennington et ah 2014)). That is, if one sorts {Xyj^yji} in 
decreasing order, then the r-th value in the list is roughly 


^[r 


k 

7*^/4 


where k is some constant. Some calculation shows that 


4fc, E 


X, 




/ <x 




and thus when a: is a small constant 

So there are only a negligible mass of Xy,^yji’s that are small constants, which vanishes when L increases. 
Furthermore, we empirically observe that the relative error of our objective is 5%, which means that the 
errors induced by Xy,^yjds that are small constants is only a small fraction of the objective. Therefore, 
Su) w' ^w,w'0{\Ayj^yji\^) is small compared to the objective and can be safely ignored. 
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